Two-channel Kondo phases and frustration-induced transitions in triple quantum dots 
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We study theoretically a ring of three quantum dots mutually coupled by antiferromagnetic ex- 
change interactions, and tunnel-coupled to two metallic leads: the simplest model in which the 
consequences of local frustration arising from internal degrees of freedom may be studied within a 
2-channel environment. Two-channel Kondo (2CK) physics is found to predominate at low-energies 
in the mirror-symmetric models considered, with a residual spin-| overscreened by coupling to both 
leads. It is however shown that two distinct 2CK phases, with different ground state parities, arise 
on tuning the interdot exchange couplings. In consequence a frustration-induced quantum phase 
transition occurs, the 2CK phases being separated by a quantum critical point for which an effective 
low-energy model is derived. Precisely at the transition, parity mixing of the quasi-degenerate local 
trimer states acts to destabilise the 2CK fixed points; and the critical fixed point is shown to consist 
of a free pseudospin together with effective 1-channel spin quenching, itself reflecting underlying 
channel-anisotropy in the inherently 2-channel system. Numerical renormalization group techniques 
and physical arguments are used to obtain a detailed understanding of the problem, including study 
of both thermodynamic and dynamical properties of the system. 

PACS numbers: 71.27.-|-a, 72.15.Qm, 73.63.Kv 



I. INTRODUCTION 

Quantum dot devices have recently been the focus of 
intense investigation, due to the impressive experimental 
control available in manipulating the microscopic inter- 
actions responsible for Kondo physics^"^. In particular 
the classic spin-i Kondo effect^ — in which a localized 
spin is fully screened by coupling to itinerant conduction 
electrons in a single attached metallic lead — has now 
been widely studied experimentally (see eg [4-8]). 

But arguably the most diverse and subtle Kondo 
physics results from the interplay between internal spin 
and orbital degrees of freedom in coupled quantum dots. 
A wide range of strongly correlated electron behavior is 
accessible in such systems, with variants of the standard 
Kondo effect^"^^, quantum phase transitions^^'^'^'^^"^^ 
and non-Fermi liquid phases^^"^^ having been studied 
theoretically. In particular, both double-^°'^^'^^'^'^"^* 
and triple^^'^^'^'^'^^"'^® quantum dots have been shown 
to demonstrate low-temperature behavior quite differ- 
ent from their single-dot counterparts^. Advances 
in nanofabrication techniques''"*''^^"^'', and atomic-scale 
manipulation using scanning tunneling microscopy 
now allow for the construction of coupled quantum dot 
structures, in which the geometry and capacitance of the 
dots can be controlled, and their couplings fine-tuned^'"^^. 
Experimental access to a rich range of Kondo and related 
physics is thus within reach. 

One of the most delicate effects, however, arises in the 
two-channel Kondo (2CK) model proposed by Nozieres 
and Blandin^'', which describes a spin-i symmetrically 
coupled to two independent metallic leads. The stan- 
dard, strong coupling Fermi liquid fixed point common 
in single- channel models is destabilized here. Much stud- 
ied theoretically (for a review see [48]), the dot spin in 



the 2CK model is overscreened at low temperatures, with 
each channel competing to compensate the local mo- 
ment. The nontrivial intermediate-coupling fixed point 
which controls the low-temperature behavior of the sys- 
tem has a number of unusual non-Fermi-liquid (NFL) 
properties^^^^^, including a residual entropy of ^ ln(2) 
and a magnetic susceptibility which diverges logarithmi- 
cally at low temperature. 

The key ingredient of this 2CK physics is the sup- 
pression of charge transfer between the two symmetri- 
cally coupled leads'''''^^'^^ — the necessity in essence of 
a central spin. Experimental realizations of 2CK physics 
have been variously sought in eg heavy fermion systems 
containing Uranium^'^^^^, scattering from 2-level systems 
using ballistic metal point contacts^^''^'', systems with 
tunneling between nonmagnetic impurities in metals^®'^^, 
and from impurities in graphene^*^; albeit that the inter- 
pretation of observed behavior in terms of 2CK physics 
is not always unambiguous. 

Recently however, Potok et al have demonstrated''''^ 
clear two-channel behavior in a coupled quantum dot de- 
vice in which one small and one large dot are tunnel- 
coupled, with the small dot also coupled to a metal- 
lic lead. The large dot acts as an interacting second 
lead, but is fine-tuned to the Coulomb blockade regime 
so that charge transfer is energetically disfavoured. A 
small degree of inter-lead charge transfer must nonethe- 
less occur, so ultimately the system is a Fermi liquid, 
with the 2CK fixed point destabilized below some low- 
temperature scale (which crossover has been widely stud- 
ied theoretically'*'^*'^^, as has the instability of the 2CK 
fixed point to channel asymmetry*^'^''^^'^^). But at finite 
temperatures, 2CK physics is clearly observed^^. 

Studying coupled quantum dot systems in a two- 
channel environment has attracted considerable theo- 
retical interest recently-' '''^^'^'''^^'^■^"^^, in part because 
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NFL states are accessible. The two-impurity two-channel 
Kondo model^^'^^'^^ (where two antiferromagnetically 
coupled dots are each coupled to their own lead) is a 
prime example, in which the tendency to form a local sin- 
glet state on the two dots competes with the formation 
of two separated single-channcil Kondo states^^'^®'^'^. The 
quantum critical point separating these phases is again 
the 2CK fixed point^*'*'''. Triple quantum dot (TQD) 
models with three dots coupled to two leads in a mirror- 
symmetric fashion have also been studied^^'^^'^*'^^'^^'^°, 
recent work^^ by Zitko and Bonca showing in particular 
that a range of fixed points familiar from simpler quan- 
tum impurity models are accessible in a ring model. In- 
deed both two-channel and two-impurity Kondo effects 
are realised, on tuning the interdot couplings as a third 
dot is coupled to the two-impurity two-channel system^^. 

TQDs arranged in a ring geometry also provide^^'^^'"*^^ 
the simplest and most direct means of studying the ef- 
fect of local frustration on Kondo physics. In mirror- 
symmetric systems, all dot states can be classified by 
their parity under left-H-right interchange^ This sym- 
metry permits a level crossing of states in the isolated 
trimer on varying the interdot coiiplings, with a pair of 
degenerate doublets comprising the ground state when 
all dots are equivalent. It was recently shown^^ that 
this level-crossing is preserved in the full many-body 
system when a single lead is attached. The situation 
is however more subtle on coupling the trimer to two 
leads, which is the focus of the present paper. We 
study a two-channel TQD ring model, shown schemat- 
ically in Fig. 1 (and discussed below), function of 
the interdot exchange couplings; using the density ma- 
trix extension'^^''^^ of Wilson's numerical renormalization 
group (NRG) technique . 

We show that 2CK physics predominates in this model, 
but that two distinct 2CK phases must in fact arise since 
local trimer states of different parity can couple to the 
leads. In consequence, one expects a quantum phase 
transition to occur between the two 2CK phases. This is 
indeed shown to arise, with the 2CK phases separated by 
a nontrivial quantum critical point, the nature of which 
is uncovered explicitly and analysed in detail. 

The paper is organized as follows. In Sec. II we discuss 
the two-channel TQD Hamiltonian, and develop low- 
energy effective models to describe the behavior of the 
system when deep in each 2CK phase. Symmetry argu- 
ments indicate the presence of a quantum phase transi- 
tion near to the point of inherent magnetic frustration 
in the TQD, and an effective low-energy model valid in 
the vicinity of the transition is derived. Sec. Ill presents 
NRG results for the full system, considering both thermo- 
dynamics and dynamical properties of the 2CK phases. 
The transition itself is investigated in Sec. IV, and the 
nature of the critical point elucidated, employing heuris- 
tic physical arguments in addition to direct calculation. 
In Sec. V the effective low-energy model describing the 
transition is itself studied directly using NRG. The paper 
concludes with a brief summary. 




FIG. 1: Schematic of the two-channel quantum dot trimer. 



Before proceeding, we point out that the stability of 
2CK physics in the model studied here is of course del- 
icate, just as it is for the standard two-channel Kondo 
model^''''^^'^^. A small degree of charge transfer, which 
would arise in a real TQD device from interdot tunnel- 
couplings (as opposed to exchange couplings), will ulti- 
mately lead to a crossover^*'^^'^^ from a 2CK to a sta- 
ble Fermi liquid fixed point, below some characteristic 
low-temperature scale. The same situation is of course 
well known to occur for channel anisotropy'*^'^^'^^'^'^, as 
would arise in the TQD model upon destruction of overall 
leftoright symmetry via eg different exchange couplings 
J23 ^ J21 (see Fig. 1). 



II. 2CK TRIMER MODEL 

Wc consider a system of three (single-level) quantum 
dots, arranged in a triangular geometry, as illustrated 
in Fig. 1. Dot '2' is coupled to dots '1' and '3', which 
are coupled to each other and to their own metallic 
lead. We model the central dot 2 strictly as a spin- 
\ to prevent inter-lead charge transfer, but dots 1 and 
3 are Anderson-like sites, permitting variable occupa- 
tion. Tunneling is allowed between these terminal dots 
and their connected leads, but the dots are coupled to 
each other by an antiferromagnetic (AF) exchange in- 
teraction to form a Heisenberg ring. We focus explic- 
itly on a system tuned to left/right mirror symmetry 
(see Fig. 1), with Hamiltonian H = Hq + i?hyb + -f^tri- 

Here Hq = Y.a,k,o- ^k^Ika^akcr ^^^^^^ ^"^^ equiva- 

lent non-interacting leads (a = L,R), which are tunnel- 
coupled to dots 1 and 3 via -ffhyb = X^k a ^i'^^Lka^la + 

^^Rka'^Sa + H.c). iftri describes the trimer itself, with 
exchange couplings J and J', 

i?tri = JS2 • (Si + S3) + j'Si • S3 

-I- e(ni +713) + U (nii-ni; -|- fis^fiai) 

where Si is a spin-i operator for dot i. For dots i = 1, 3, 

= J2a "i<T = 4cr^ia number operator, e is 

the level energy and U is the intradot Coulomb repul- 
sion. The full Hamiltonian H is thus invariant under 
simultaneous 1 O 3 and L ^ R permutation. 
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The isolated trimer 



B. Effective low-energy models 



We are interested in the TQD deep in the Af = 3 
Coulomb blockade valley, where each dot is in prac- 
tice singly occupied. To this end we consider explicitly 
e = —^U, the A/" = 2 or 4 states being much higher 
(~ ^U) in energy. iJtri (and the full H) is then particle- 
hole symmetric, although this is incidental: we require 
only that the singly-occupied manifold of TQD states 
lies lowest. This manifold comprises two doublets and a 
spin quartet, which is always max( J, J') higher than the 
doublets for AF exchange couplings. 

For any J, J', the lowest doublets of the isolated TQD 
are 
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with = I and cr = ± for spins t/4-, and = 
the raising/lowering operator for the spin on dot 2. 
|vac) = J2a2 l~'^2; — ) defines the 'vacuum' state of the 
local (dot) Hilbert space, in which dots 1 and 3 are un- 
occupied, and dot 2 carries a spin-i. 

The energy separation of the two doublets is E'a = 
— = J — J' , with |-|-; S") the ground state of the 
isolated TQD for J' > J and |-; 5^) lowest for J' < J. 
When J' is dominant, dots 1 and 3 naturally lock up into 
a singlet (see l-l-;^^), Eq. 2), leaving a free spin on dot 
2; with (Si . S3) = -| and (Si • S2) = (S3 • 83) - 0. 
For J' < J by contrast, dots 1 and 3 are now in a triplet 
configuration (|-; S^) in Eq. 2), with (Si • S3) = +j and 

(Si • S2) = (S3 • S2) = -i 

The states are degenerate precisely at J' = J, reflect- 
ing the inherent magnetic frustration at that point. A 
level crossing of the doublets at J' = J is permitted be- 
cause each has different symmetry under 1 -H- 3 permu- 
tation. We define a parity operator Pij which exchanges 
orbital labels i and j (as discussed further in the Ap- 
pendix). From Eq. 1 it is clear that Pi, 3 commutes with 
the isolated TQD Hamiltonian, [Pi,3,^?tri] = 0. Thus all 
states of Htri can be classified according to 1 O 3 parity, 
the eigenvalues of Pi. 3 being ±1 only (since Pf 3 = 1). In 
the spin-only ('singly-occupied') regime, the parity oper- 

ator may be expressed concisely as^^ Pi^s = (1 — S13), 
with S13 = (Si -I- S3) the total spin of dots 1 and 3. Thus 
Pi^3\±;S^) = ±\±;S^) describes the parity of the dou- 
blet states of i/tri- The full lead-coupled Hamiltonian H 
is not of course invariant to 1 <s-> 3 interchange alone, but 
rather to simultaneous exchange of the 1, 3 dot labels and 
the left/right leads (embodied in Pl,r = Pi, 3 Ilk ^ik,iikj 
see Appendix); which we refer to as 'overall L R" sym- 
metry. 



On tunnel-coupling to the leads, effective models de- 
scribing the system on low-cnergy/temperaturc scales in 
the N = 3-electron valley of interest may be obtained 
by standard Schrieffer- Wolff transformations^'''^, pertur- 
batively eliminating virtual excitations into the Af = 2- 
and 4-electron sectors of H to second order in i/hyb (and 
neglecting retardation effects^ as usual). The calcula- 
tions are lengthy, so rather than giving full details we 
sketch below a somewhat simplified, but physically more 
transparent, account of the key results (Eqs. 7,13 below). 

First we consider the effective low-energy model ap- 
propriate to the temperature range {J, J') < T < U, 
in which all dots become singly-occupied. Here the ap- 
propriate unity operator for the TQD Hilbert space is 
iso = E,Ti,<T2,<T3 |o-i;o-2;o-3)(o-i;f72;f73|, with ai =t/4- the 
spin of dot i. To second order in the dot- lead tunneling 
V, a spin-model of form Hg = Hq + H^g + H^^.^ arises; 
where Hq describes the leads as above and: 



= JS2 • (Si + S3) + J'Si • S3 



(3) 



Here the effective exchange coupling is pJx = 8T/{'kU), 
with p the lead density of states per orbital at the Fermi 
level; and F = vrV^p^ is the hybridization, with total 
lead density of states pj, = Np, and A'' the number of 
orbitals/k-states in the lead (such that NV^, and hence 
r, is finite in the continuum limit TV — > 00). In Eq. 3, s^q 
is the spin density of lead a = L{R) at dot 1(3) given by 



:Kk(T ' 



(4a) 
(4b) 



with cr the Pauli matrices and fJ^Q^ the creation operator 
for the 'O'-orbital of the a = L,R Wilson chain. 

As above, the lowest states of H^^.^ are the doublets 
|±;S'^) given in Eq. 2. Provided they are not near- 
degenerate, only the lower doublet need be retained: 
\+; S*^) for J' > J and |-; S^) for J' <C J. To first order 
in Jk, an effective low-energy model is then obtained 
simply by projecting into the reduced Hilbert space of 
the lowest doublet, using 



i±=^|±;5^)(±;5^ 



(5) 



for the appropriate 7 = -|- or — doublet ground state. 
The resultant Hamiltonian i?eff,7 = ^-yH^f^l^ follows as 

HeS,-y = yxi-yiSl + 83)1^ • (SLO + Srq) (6) 



using the symmetry l-ySil-., 



1^831^. 
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Eq. 6 is of two-channel Kondo form, 



i-ySii-y' = 77'i-yS3i-y', 0116 obtams 



eff,7 



= J. 



S7 • {slo + sro), 



(7) 



where is the spin-i operator representing the ap- 
propriate doublet 7 = + or — , with components = 

\T,S')S^r,S'\ and S± = I7; ± i)(7; t^I- At this 
level of calculation the effective exchange coupling is 
given by Jkj = (7; + S^lr, +h)JK; so from Eq. 2 

an AF effective Kondo coupling Jk- = +^Jk then arises 
for the | — ; S^) ground state appropriate to J' <^ J, while 
for the |+; 5^) doublet (lowest for J' > J), Jk+ = re- 
sults. In the latter case, there is in fact a weak residual 
AF coupling: a full 0{V^) Schrieffer- Wolff calculation 
gives precisely the two-channel Kondo model Eq. 7 as 
one would expect, but with Jk± given by 



pJk+ 



_4T r 1 

~V 1 3( J' - J 



1 



pJk- =- 



4r 



J) + 2U 3J' + J + 2U 

9 1 



377 {{J- J') + 2U 5J-J' + 2U 



(8a) 

; (8b) 



yielding pJx- = ^pJk oc ^ to leading order in and 
a much smaller but non- vanishing pJk+ — + C'(f73') 
for the singlet-locked doublet |-|-; S^), reflecting the resid- 
ual AF coupling between the spin on dot 2 and the leads. 

Hence, sufficiently deep in either regime J' ^ J or 
J' ^ J, the low-energy behavior of the system is that of 
a 2CK model. The lowest spin-i state of the TQD is thus 
overscreened by conduction electrons, embodied in the 
infrared 2CK fixed point describing the non-Fermi liquid 
ground state^^'"'*'^^, in which the partially quenched spin 
is characterised by a residual entropy of Simp = ^ In 2 
{ks = 1); overscreening setting in below the character- 
istic two channel Kondo scale , determined from 
perturbative scaling as^'' 



^Ivf ~ Jif^exp(-l/pJK-y). 



(9) 



Since Pi 3 commutes with all components of S-^ in Eq. 7, 

[^fcfi,7, A,3] = 0, whence 1 O 3 parity is conserved in 
the effective low-energy model. Since that parity is deter- 
mined by 7 = -|- or — , there are two distinct 2CK phases, 
which one thus expects to be separated by a quantum 
phase transition (QPT). 

In the vicinity of the transition, ie close to J' ~ J, nei- 
ther of the two 2CK models in Eq. 7 is of course sufficient 
to describe the low-energy physics: the states \+;S^) 
and | — ; S^') are now ncar-dcgencratc, so both must be re- 
tained in the low-energy trimer manifold. Hence, defining 
1 = 1+ -h 1- and proceeding in direct parallel to the dis- 
cussion above, an effective low-energy model in the vicin- 
ity of the transition is obtained from = l(-ffcff + 
H^,i)^. From Eq. 3 for H^^^, using Si = A,3S3A,3 such 
that (7;S'^|Si|7';S'^') = 77'(7; 5^|S3|7'; -S^') and hence 



^^eff"' = E 5"^^[i7(Sl + S3)i-y " (§L0 + Sro) 

+i^(Si - Sa)!--, • {SLO - SRo)] (10) 

+i£A (i+ - i-) . 

The final term in Eq. 10 (arising from IH^^^l) is simply 
the energy difference between the 7 — +/— doublets, 
and the first term is (see Eq. 6) the 2CK coupling of each 
doublet to the leads. It is helpful to recast Eq. 10 in 
terms of spin-i operators for real spin S (= S-y) and 

pseudospin 7" for the local Hilbert space, defined by 



and 



7 



(11a) 
(lib) 

(12a) 
(12b) 



From Eq. 11, the TQD doublets I7; 5^) arc each cigen- 
states of and T^; in particular, the eigenvalues of 
= |(1+ — 1-) correspond simply to (half) the par- 
ity of the appropriate doublet. By contrast, the doublets 
are interconverted by 7"=^ (Eq. 12b), T+|-; 5^) = \+; S") 
and T~\+; S^) = \ — ;S^), acting to switch parity. 
After simple if laborious algebra, Eq. 10 reduces to 

Hlr' = (I Ja + Jst^)S • {SLO + SRo) 

+ JmUf+ + f-)S ■ {SLO - Sro) (13) 

in terms of spin/pseudospin operators, with 

Ja = {Jk+ + Jk-) Jb = {Jk+ - Jk-) (14) 

and Jmix = {+]+\\Sl - S^\-]+\)Jk = -^Jk using 
Eq. 2. [A full 0(y2) Schrieffer- Wolff calculation again 
gives precisely the effective Hamiltonian Eq. 13, with 



pJn 



2r 



+ 



VStt 1 3 J' + J -(- 2f/ U-J' + 2U 



+ 



+ 



3iJ'-J)+2U {J-J') + 2U 



(15) 



recovering Jmix = '^Jk to leading order in l/U, and 
Jk± given by Eq. 8.] 

Eq. 13 is the essential low-energy model applicable 
to the vicinity of the QPT; we study it directly via 
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FIG. 2: Lowest five NRG energy levels of the left-charge/right-charge/spin subspaces [Ql, Qr, S") = (0, 0, i) [solid lines] and 
(0, 1, 0)/(l, 0, 0) [dashed], versus (even) iteration number, A''. Shown for pj = 0.1, U = U/{nT) — 7, with pj' = [panel (a)] 
and pJ' = 0.105 [(b)]. For comparison, (c) shows corresponding results for a pure (single-spin) 2CK model with pJx = 0.055. 



NRG in Scc.V. The pscudospin operators can natu- 
rally be classified according to 1 3 parity, the com- 
ponents of T having different parity under Pi, 3: using 
Pi,3\t,S'){7';S^\ = 77'l7;S^>(7';S1Pi,3, Eqs.^ llb,12b 

give A.st^ = +t^A,3, while A,3t± = -r±A,3 for 
the raising/lovi^ering components. By contrast, all com- 
ponents of spin S (Eqs. 11a, 12a) commute with A,3- 
Hence, since the global (overall L ^ R) parity must be 
conserved ([il*^"^, A,fl] = 0), interactions involving the 
component of pseudospin can only couple to even 
combinations slo + of the lead spin densities (sym- 
metric imdcr interchange of the L/R lead labels), as in 
the first term of Eq. 13 (or 10); while by the same to- 
ken interactions involving T"^ must be associated with 
the odd (antisymmetric) combination s^o — ^ in 
the second term of Eq. 13. In the vicinity of the QPT 
the latter is of course the key interaction in iJ^g-"^, since 
in switching the parity of the TQD states it in essence 
drives the transition between the two 2CK phases. 

Finally, note that the last term of Eq. 13, equivalent to 

a magnetic field acting on the pseudospin, energetically 
favors the \ — ;S^) doublet (T^ = ~\) when their energy 
separation Ea > 0, and |-h;S"=) (T'^ = -t-^) for Ea < 0. 
Hence, when \Ea\ is sufficiently large that only one of 
the doublets need be retained in the low-energy TQD 
manifold, the terms are obviously suppressed; and 
Eq. 13 then reduces as it must to one or other of the 
2CK models Eq. 7. In fact, as shown in Sec. IV ff, for 
any |i<^A| 7^ 0, the EaT^ term in Eq. 13 ensures that one 
or other of the 2CK fixed points ultimately remains the 
stable low-temperature FP. 

By contrast, precisely at the point of frustration where 
|£'a| = 0, the terms in Eq. 13 destabilize the 2CK 
FPs. The resultant H^^^^ then describes the quantum 
critical point which separates the two 2CK phases, and 
which we consider in detail in Sec. IV. 



III. PROPERTIES OF THE 2CK PHASES 

The physical picture thus indicates that 2CK physics 
dominates the low-energy behavior of the model; with 
a QPT occurring as a function of J' between two 2CK 
phases of distinct parity. 

We now analyse the properties of each 2CK phase of 
the full model; using Wilson's NRG technique^^'''*, em- 
ploying a complete basis set of the Wilson chain^'^ to 
calculate the full density matrix^^'^^ (for a recent review 
see [75]). Calculations are typically performed for an 
NRG discretization parameter A = 3, retaining the low- 
est Ns = 4000 states per iteration. As above we choose 
for convenience e = —\U, and consider a symmetric 
constant density of states for each lead, with density of 
states per conduction orbital p = 1/(2D), and bandwidth 
D/T = 100 » 1 (such that results shown are essentially 
independent of D). In all calculations shown explicitly, 
we use a fixed pJ = 0.1 and U = [//(ttF) = 7, varying 
the exchange J' (see Fig. 1)''^. 

Fig. 2 shows the evolution of the lowest energy lev- 
els of the system as a function of NRG iteration number 
TV, exemplifying RG flow between different FPs of the 
model. Panel (a) is for a system deep in the J' < J 
regime (specifically pJ' = 0), while panel (b) shows the 
energy levels for pJ' = 0.105 (> pJ). For comparison, 
panel (c) is for a pure (single spin-^) 2CK model of form 

Eq. 7, ie H = J^S ■ (s^o + §^0)7 with a Kondo couphng 
pJk = 0.055 chosen to be the same as the effective cou- 
pling pJk- of the ground state TQD doublet in (a) (as 
obtained from Eq. 8(b)). In both cases (a) and (b), the 
levels are seen to converge quite rapidly to their iV — > 00 
values, which are clearly those of the 2CK FP in Fig. 2(c). 
These levels are of course characteristic of the 2CK FP, 
and — after a trivial rescaling by a factor which depends 
on the NRG discretization parameter A arc described 
by the fractions 0, ^, ^, |, 1... as determined from a 
conformal field theory analysis of the FP^^'^^. 

The iteration number by which the levels converge to 
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FIG. 3: Thermodynamics of the J' < J phase. Shown for 
fixed U — 7 and pj — 0.1, varying pj' = (sofid lines), 0.0375 
(dotted) and 0.075 (dashed). Panel (a) shows 5'imp(T), (b) 
shows Tximp(T), both vs T/F. Insets show the universality 
arising on rescaling in terms of the 2CK Kondo temperature, 
Tk, all curves collapsing to a common form. 

the set of 2CK FP levels is however clearly different in (a) 

and (b) , reflecting the different Kondo scales characteris- 
tic of the two cases. Case (b) (J' > J) flows close to the 
local moment (LM) FP between .A^ = 10 to 60 (with lev- 
els naturally characteristic of the LM FP in that range), 
and approaches the stable 2CK FP by -ZV ~ 80. By con- 
trast, convergence to the 2CK FP in (a) and (c) both 
occur at -/V « 50 (with a much shorter range of N close 
to the LM FP). Since the iteration number is related to 
an effective temperature through^''''^'* T/F ^ A~^/^, the 
2CK Kondo scales of the three examples in Fig. 2 are 
thus exponentially small (as expected from Eq. 9); but 
with /T for (b) being some 8 orders of magnitude 

smaller than that of (a), reflecting the distinct nature of 
the coupling in the two cases, as discussed in Sec. II B. 

A. Thermodynamics 

NRG results for thermodynamics in each phase are 
now considered. We focus on the 'impurity' (TQD) 
contribution^^'^^ to the entropy, S'imp(r), and the uni- 
form spin susceptibility, Ximp(r) = {{S^)^)inip/T (here 
refers to the spin of the entire system and (n)iinp = 
(ii) — (f2)o with (0)o denoting a thermal average in the 
absence of the TQD); the temperature (T) dependences 




10"* 10"" 
T/F 



FIG. 4: Spin susceptibility Ximp(T) ?;s T/T in the J' < J 
regime, for same parameters as Fig. 3. Inset: scaling collapse 
onto the universal form TK-Ximp(T') = aln(T/f/T). 



of which provide clear signatures of the underlying FPs 

reached under renormalization on progressive reduction 
of the temperature/energy scale''^'^^. 



1. J' <J 

For the J' <^ J regime the effective 2CK model Eq. 7 
should describe the system for T <^ \J — J'\ = I-EaIj 
where the lowest TQD doublet, in this case the odd parity 
state | — ;S'^), couples symmetrically to the leads. 2CK 
physics is thus expected below T ~ Tk = T]fJ^ . 

For the full model. Fig. 3 shows both S'imp(T) [panel 
(a)] and Tximp{T) [panel (b)] for fixed [/ = 7 and 
pJ = 0.1, with J'/ J = 0, 0.375 and 0.75, aU deep in 
the J' < J regime. In each case, at high temperatures 
T > U the behavior is governed by the free orbital^"* 
(FO) FP, with all possible states of the TQD thermally 
accessible, and hence Simp = ln(32). For T %U the dots 
become in essence singly occupied, and hence an entropy 
of ln(8) is expected^°. Below T ~ |.Ea| all but the lowest 
trimer doublet is projected out. Thus S\mp approaches 
ln(2), signifying the LM FP where the lowest doublet 
is essentially a free spin-i disconnected from the leads 
(Eq. 7 with Jk- = 0). 

The LM FP is however unstable, and the system al- 
ways flows to the stable 2CK FP below T ^ Tk, recov- 
ering (Fig. 3(a)) the | ln2 residual entropy known to be 
characteristic of the 2CK PP'^^"^^. In practice we may 
deflne a Kondo temperature Tk = T^ through the en- 
tropy, via S'impCT^) = 0.55 (suitably between In 2 and 
i In 2) : or alternatively Tk = T^ through the spin sus- 
ceptibility via T^XimpiT]^) = 0.07 (as in [74]). Deep in 
the 2CK phases the two definitions arc of course equiva- 
lent (T^ = = Tk), probing as they do the common 
characteristic scale associated with flow to the 2CK FP. 
The inset to Fig. 3(a) shows the entropy of the three sys- 
tems rescaled in terms of T/Tk; showing scaling collapse 
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FIG. 5: Thermodynamics of the J' > J phase. Shown for 
fixed U = 7, pJ = 0.1, varying pj' = 0.105 (solid lines), 0.115 
(dotted) and 0.125 (dashed). Panel (a) again shows S'imp(r), 
(b) shows rximp(T), both vs T/T. Insets show rescaling in 
terms of the Kondo temperature, Tk, yielding the same uni- 
versal 2CK scaling curves as in Fig. 3. 



to a common functional form, ie the universality char- 
acteristic of the crossover from the LM FP to the stable 
2CK FP^K 

The underlying FPs of the model are likewise evi- 
dent from Tximp(T), Fig. 3(b). The highest T behavior 
corresponds to two uncorrclatcd sites (dots 1 and 3 of 
the TQD) and a free spin (dot 2). Hence TximpiT) = 
2 X I + I = i, readily understood as the mean {{S^)^) 
of the quasidegenerate states. On decreasing T the LM 
FP is again rapidly reached, the lowest TQD doublet fol- 
lowing the free spin-^ Curie law TximpiT) = j. Below 
Tk the spin susceptibility is quenched'*^'^-'^'^'^ in the sense 
that Tximp(r) -)> as T ^ 0. In the inset to Fig. 3(b) 
the data are rescaled in terms of T/Tk, again showing 
universality in the approach to the 2CK FP^^. 

The characteristic low-T logarithmic divergence^®'^^'^^ 
of Ximp(T') itself is evident in Fig. 4 for the same param- 
eters as Fig. 3. The slopes of the divergence vary widely 
between the three cases, but all collapse to the universal 
form'"'^ TKXimp{T) = a\n{TK /T) (with a a constant), as 
seen in the inset of the figure. 



2. J' > J 

Prom the arguments of Sec. II B, 2CK physics is again 

expected at low-energies in the J' ^ J regime, where the 
even-parity |+; S"^) doublet is now the lowest TQD state. 
This is confirmed in Fig. 5, where (in direct parallel to 
Fig. 3 for J' <^ J) thermodynamics are shown for fixed 
U = 1 &ud pJ = 0.1, but now with J' / J = 1.05, 1.15 
and 1.25: the scaling forms for Simp(T') and Tximp{T) vs 
T/Tk, shown in the insets of Fig. 5, are precisely those 
arising for J' <^ J and shown in Fig. 3. 

While the low-energy physics in the two regimes J' ^ 
J and J' ^ J have common 2CK scaling behavior, we 
add that the Kondo scales Tk themselves evolve very 
differently with J'; as will be considered explicitly in 
Sec. IIIC (see Fig. 10), and evident even from cursory 
comparison of Figs. 5,3. 



B. Dynamics 

We turn now to dynamics, focussing largely on the 
local single-particle spectrum of dot 1 (or, equivalently 
by symmetry, dot 3): Di(oj) = —^lmGi{uj), with Gi{uj) 
the local retarded Green function. We obtain it through 
the Dyson equation, 



[Gi(o 



-1 



[G;(a;)]-^-Si(a;), 



(16) 



where G''1{lu) is the non-interacting propagator (ie for 
U = = J = J'); and Si(a;) is the proper electron 



self-energy, with T,i{uj) = (w) — iT,{{ui) (such that 
(w) > 0). The non- interacting G5(a;) is given trivially 
byi [G?(cj)]-i = UJ+ -€- r(w) (with uj+ = uj + iO+), 
where T{uj) = J2\i ^^[w"'" — ek]~^ is the usual one-electron 
hybridization function for coupling of dot 1 to the left 
lead: r(w) = r^{uj)-ir^ {u), with r^(w) = F (= nV^Pj,) 
for all |a;| < D inside the band, and T^{uj = 0) = 0. 

An expression for the electron self-energy Ei (cj) is also 
readily obtained using equation of motion methods^'^^. 
It is given by 

Ei(u;) = [Gi(a;)]-i|c/((c,^ni^;cl^)), 



(17) 



where {{A; B))^ is the Fourier transform of the retarded 

correlator {{Aih); B{t2))) = -ie{h-t2){{Aiti), B{h)}), 
and where the local Green function itself is given by 
Gi(a;) = {{ci^', Ci^)) u (independent of spin a in the ab- 
sence of a magnetic field). The self-energy can be cal- 
culated directly within the density matrix formulation 
of the NRG^i^^2,75,82^ 17; with Gi(w) then ob- 

tained from Eq. 16. Di{oj) calculated in this way is highly 
accurate^^'^^, and automatically guarantees correct nor- 
malization of the spectrum^^'^^. 
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To motivate study of the spectrum Di (= Di{u};T)), 
we note that it controls the zero-bias conductance 
through dot 1. The TQD does not of course mediate cur- 
rent from the L to R leads, since the internal couplings 
between constituent dots are pure exchange (Fig. 1 and 
Eq. 1). However the L and R leads in Fig. 1 can obviously 
each be 'split' in two (symmetrically, to preserve over- 
all i O -R symmetry), enabling a current to be driven 
through dot 1 or dot 3; with the same zero-bias conduc- 
tance in either case, by symmetry. Following Meir and 
Wingreen^^, the resultant conductance follows as 



Go 



dfjco) 



7TrDi{uj;T) duj (18) 



where f{uj] 



1 



is the Fermi function (with 



CO = the Fermi level) and Go = 2e^ /h the conductance 
quantum. The (w, T)-dependence of Di thus controls 
the conductance, and for T = in particular Gc{T = 
0)/Go = nTDiiuj = 0;T = 0). From Eq. 16, the 
local propagator for w = = T may be expressed as 
[Gi{u) = 0)]~^ = — e* -|- ir* in terms of the renormalized 
single-particle level e* and renormalized hybridization F*, 
given by 



e* = e + Sf(0) 

F* = r-t-s((o) 



(19a) 
(19b) 



in terms of the lj = value of the self-energy at T = 0; 
and hence from Eq. 18: 



GcjT = 0) F 1 

= 7rrDi(0;0) = — ——^ 

Go J- l+(fqr) 



(20) 



NRG results for single-particle dynamics are consid- 
ered below, but first a question arises. We have argued 
above that, sufficiently deep in either regime J' < J or 
J' > J, the low-energy physics of the full three-site TQD 
model must reduce to that of a single spin-i 2CK model 
(of form Eq. 7). The question is; to which dynamical 
property of a pure 2CK model should the spectral den- 
sity '}tTDi{uj) be compared? To answer this, note first 
that Eq. 18 may be written equivalently as 



GcjT) 
Go 



f 

J — C 



dm 



-Trpjnim iL(w) dw (21) 



in terms of the t-matrix, tL^co), for the L lead; with ^0(0;) 
defined in the usual way in terms of scattering of electron 
states in the a = L,R lead, via 



G. 



ak.ak' 



(a;) = 



Ck 



+ 



1 



Ck 



tocioj) 



Ek' 



(22) 



where Gak,ak'{uj) = ((Oak^; «akv))" is tii° propagator 
for the lead states. Using equation of motion methods"'^'^"'^ 
it is straightforward to show that tLiuj) = V'^Gi{ljj) (like- 
wise i/{(w) = V^G3{oj) = tL{ijj)); so —Trprplm.tL{oj) = 



7rFDi(w;T) (recall F = nV'^pj,), hence the equivalence 
of Eqs. 18,21. 

To compare 'kTDi{uj;T) for the full TQD model to 
results for a single-spin 2CK model^^'^^'*'*"^'^, we thus 
require tL(w) = tjf'^{uj) for the latter. Using the def- 
inition of the 0-orbital of a Wilson chain (Eq. 4b), the 
Hamiltonian for a single-spin 2CK model is 



^ako- 



H2CK = X] Ekolka^o 

+ E W [^"^(4kT«ok'T - «lk;a„k';) (23) 



akk' 



and for this case equations of motion again yield Eq. 22 
for Gak.akK"^)! ^ut now with tL{u)) = -^Gs{oj) where 

G,(a;) = {{S-hoi+S'h,^;S+fl^+S^fl^))^. (24) 

Comparison of —T:prplmtL{uj) for the full TQD model 
with its pure 2CK counterpart (using Pj,/N = p = 
1/(2D)) then gives the desired correspondence 

'kTDi{uj) <H- D Pj^{uj) (25) 
with spectral density 

p^(w) = -l{pJKfluvGs{uj) (26) 
for the single-spin 2CK model**. 

1. Dynamics: results 

Fig. 6 shows the T = spectrum 'kTDi{uj) vs w/F, in 
the J' < J regime (for the same bare parameters as Fig. 3 
for thermodynamics). The main panel shows results on 
a log scale for u; > 0; and particle-hole symmetry for 
e = — if/ means Di{(jj) = Di{—ui) (as seen in the inset). 
The figure also shows Dpk{oj) for the single-spin 2CK 
model (for pJk = 0.07, ensuring an exponentially small 
Kondo scale but otherwise chosen arbitrarily). 

Wc first comment on the high-energy spectral features 
('Hubbard satellites') in Di{u). As usual, these reflect 
simple one-electron addition to the isolated TQD ground 
state. The lowest such excitation from the | — ;S'^) TQD 
ground state incurs an energy cost AE- = j{J — J') — e. 
On tunnel coupling to the leads these features are nat- 
urally broadened, but the satellites are centered on a; ~ 
AE-, indeed seen from Fig. 6 to shift slightly to lower 
uj on increasing J'. In the J' > J regime (not shown), 
directly analogous behavior arises. Here the \+;S^) dou- 
blet is the TQD ground state, the Hubbard satellites in 
Di{ijj) are now centered around cj ~ AE^ ~ |(J' — J) — e, 
and hence shift to higher frequency as J' is increased. By 
contrast, see Fig. 6, Hubbard satellites are simply absent 
in Pi^{lo) - addition/removal excitations are suppressed 
by construction in modeling a dot strictly as a spin. 
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FIG. 6: T=0 local spectrum 7rrDi(cj) vs w/F in the J' < J 
regime. Shown (as in Fig. 3) for fixed U = 1 and pJ = 
0.1, varying J' /J = (solid line), 0.375 (dotted) and 0.75 
(dashed). The dot-dashed line is the spectral density Dpi^ioj) 
(Eq. 26) for a single-spin 2CK model (Eq. 23), with pJx = 
0.07. Inset: results shown on a linear frequency scale. 

The most important characteristic of the spectra in 
Fig. 6 is of course the low-energy Kendo resonance, the 
form of which reflects RG flow in the vicinity of the stable 
2CK flxed point. At the Fermi level cj = in particular, 
'kTDi{lj = 0) = ^ in all cases (likewise Dp^{0) — i) 
- ie reaches h,alf the unitarity limit, a hallmark of the 
2CK FP, likewise known from study of the two-channel 
(quadnipolar) single- impurity Anderson model®^'^^. 

Universal scaling of single-particle dynamics is con- 
sidered in Fig. 7, where '!tTDi{uj) vs uj/Tk is shown 
(with Tk defined in practice from thermodynamics, as 
in Sec. Ill A). The figure includes the three examples 
shown in Fig. 6 for the J' < J (odd parity ground state) 
regime, as well as spectra for the J' > J (even parity) 
regime for the same parameters as Fig. 5; together with 
Dpk{uj) vs uj/Tk for the pure 2CK model. As seen 
from Fig. 7, all spectra collapse perfectly to a universal 
scaling form; confirming that the TQD model - be it in 
the J' > J OT J' < J regime - is described by the 2CK 
model at low-energies. 

As shown in Fig. 7, the leading low-frequency asymp- 
totics of the scaling spectrum are found to be 

M 

nTD^iu) ^[l-a{\u\/TK)h (27) 

with the constant a ~ 0.47; as consistent with that found 
for the quadrupolar Anderson impurity model^^, and in 
contrast to the [1 — a(a;/Tff )^] decay characteristic of the 
normal Fermi liquid FP in single-channel models. For 
\oj\/Tk ^ 1 by contrast, as also shown in Fig. 7, the 
asymptotic behavior is 

M 

7rr£>i(a;) ^"r^"" A/[lnWcu\/TK) + B] (28) 

(with A and B pure constants 0(1)); which form is also 
asymptotically common to the single-channel spin-i An- 




lj/Tk 

FIG. 7: T = scaling spectrum for the 2CK phases. 7rrDi(a)) 
vs w/Tk is shown both in the J' < J regime (for the bare 
parameter sets used in Fig. 3) and for the J' > J regime 
(with the parameters used in Fig. 5). Dpxii^) vs uj/Tk for 
the single-spin 2CK model is also shown. All spectra collapse 
to the same universal form. The dotted line shows the low- 
\uj\/Tk asymptotic behavior irTDi{u;) = i[l — a(|a;|/rif ) 2], 
while the dashed line describes the high-|w|/rK- scaling be- 
havior, 7i-rDi(w) = A/[\n'^{\u}\/TK) + B]. 



derson model*'' and its SU(2N) generalization^^. This is 
physically natural, the origin of the 'high' energy leading 

^ 1/ ln^(|a;|/T/f) logarithms being spin-flip scattering-'^, 
which occurs for both single- and two-channel models. 

We consider now the electron self-energy of dot 1 in the 
full TQD model, the real and imaginary parts of which 
are shown in Fig. 8 and naturally exhibit scaling in terms 
of uj/Tk- For S{(a;) we find the Iow-uj/Tk asymptotic 

behavior to be S{(w)/r = 1 + c|w/T|P^|2 as shown 
explicitly in the insert to panel (a). At the Fermi level 
in particular, S{(a; = 0) = F, in marked contrast to a 
Fermi liquid phase for which = 0) = generically; 

and leading to a renormalized hybridization (Eq. 19b) 
r* = r + i;{(u; = 0) = 2r. The real part of the self- 
energy is shown in panel (b), and at the Fermi level in 
particular is found to be Sf (w = 0) = — e, such that the 
renormalized level e* = (Eq. 19a). 

Since F* = 2r and e* = 0, it foUows from Eq. 20 
that the T = zero-bias conductance for the TQD 
model reduces to one-half the unitarity limit, ie Gc{T = 
0)/Go = ^, as known to be the case for a single-spin 2CK 
model^^''"'^'*'*'*^. Moreover, although the explicit case for 
which we have shown results is e = —U/2 where the TQD 
model is particle-hole symmetric, we note that F* = 2r 
and e* = - and hence 7rrDi(a;) = ^ and Gc{0)/Go — ^ 
- is also found to hold robustly on moving away from 
particle-hole symmetry. 

The results above refer to dynamics for T = 0, and we 

now touch on T > 0. Universality implies 7rrDi(a;;T) 
depends on uj/Tk and T /Tk- Without loss of generality 
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FIG. 8: Electron self-energy Ei(w)/r vs ui/V for the same 
parameters as Fig. 6. Panel (a): imaginary part, E{(a;)/r; 
panel (b); real part, Ef(a;)/r. Insets: universal scaling forms 
as a function of ijj/Tk-, with the dot-dashed lines showing the 

low- frequency asymptotic behavior S{(a;)/r = c\i^ /Tk\^ , 

and Ef(w)/r = -t/T - sgn{w)c'\t^ /Tk\'^ . 

it may be cast as 

(29) 

where the scaUng function g satisfies g{Q; 0) = 1 (such 
that Eq. 29 for T = reduces asymptotically to Eq. 27 for 
7rr£'i(cLi; T = 0)). For the Fermi level w = in particular, 
we find from NRG that the leading \ow-T/Tk dependence 
of 7rr£)i(0;r) is 

T 

7rr£»i(0;T) i[l - 6 (T/Tk)5] (30) 

with 6 ~ a a constant. Employing Eq. 29 in Eq. 18 
for Gc{T), and rescaling the resultant non-trivial integral 
(arising from the second term of Eq. 29) by exploiting the 

fact that the Fermi function depends solely on cj/T, gives 
the leading Iow-T/Tk behavior of the zero-bias conduc- 
tance as 

^ = i[l- 7 (T/T^)^] (31) 

where 7 is a constant {j = b + a' with a' = a^/irrji^) ~ 
1.07a and r]{z) the Dirichlet ry- function) . And the be- 
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uj/T 



FIG. 9: Local dynamic susceptibility FxiC^^) vs uj /F at T = 
for the same bare parameters as Fig. 6. Inset: showing 
universal scaling of Tr:Xi(<^) i^/Tk- 



havior Eq. 31 is precisely that known to arise for the 
single-spin 2CK model^^-^^.s^.sT 

Finally, we consider briefly the local dynamic spin 
susceptibility of dot 1 in the 2CK phases, given by 
Xi(u;) = -ilm((5f;5f))^. Fig. 9 shows Xi(w) vs u/T 
(for T = 0) in the J' < J regime, with the same bare 
parameter sets as Fig. 6 (the behavior discussed below 
being applicable in both 2CK phases). In the standard 
single-channel Anderson model, Xii^) exhibits character- 
istic \ow-Lu Fermi liquid behavior^'^^, Xii^) ^ ^- By con- 
trast, it is known®'5'^2 that at the 2CK FP, Xi(0) plateaus 
at a finite constant; itself proportional to the slope of the 
log divergence of XimpiT) as T ^ 0. In Fig. 4 (inset) 
for the full TQD model we showed the uniform static 
spin susceptibility to have a low-T log divergence of form 
Ximp{T) = ^ IjiItk/T), with slope a/Tx- Thus we also 
expect TkXiI^) to exhibit universality as a function of 
co/Tk; as indeed confirmed in the inset of Fig. 9. 



C. Kondo Scales 

The energy scale on which 2CK physics emerges in 
the full TQD model is of course the 2-channel Kondo 
temperature. Tk itself varies markedly with the bare 
model parameters; with qualitatively different behavior 
for J'/ J < 1 and J' /J > 1 that reflects the distinct 
nature (Sec. II B) of the lowest-energy TQD doublets in 
the two regimes, \ — ;S^) and \+;S^) respectively. This is 
shown in Fig. 10, considering explicitly the Kondo scale 
Tk = obtained (as in Sec. Ill A) from the uniform 
susceptibility, and showing T^/T vs J'/ J. In the J' < J 
regime a relatively modest increase in Tk occurs on in- 
creasing J' /J; while for J' > J by contrast Tk dimin- 
ishes very rapidly, reflecting physically the weak residual 
AF coupling (Sec. II B) between the spin on dot 2 and 
the leads, when the singlet-locked |-|-; S^) doublet is the 
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J'/J 

FIG. 10: Kondo temperature T^jV (points, with solid line 
as guide to eye) va J' / J for fixed U — 7 and pj — 0.1; as 
defined (Sec. Ill A) via T^Xin,p{T^) = 0.07. Kondo scales 
obtained from perturbative scaling (Eq. 9) are also shown, 
with pJk+ (for J'/J > 1) from Eq. 8(a) (dashed line) and 
pJk- (for J'/J < 1) from Eq. 8(b) (dotted); with a trivial 
constant prefactor adjusted to fit the data. 

TQD ground state. 

In Sec. II B wc showed that the full TQD H maps 
onto an effective single-spin 2CK model (Eq. 7) in each 
of the J' < J and > J regimes, provided the separa- 
tion between the TQD doublets \Ea\ = \ J — J'\ is suffi- 
ciently large that one or other alone dominates the low- 
energy physics. The resultant Kondo temperatures for 
the two regimes are then given from perturbative scaling 
by Eqs. 9,8. These are compared directly to the NRG 
results for in Fig. 10, and are seen to agree quantita- 
tively for \Ea\/J ^ 0.1, ie throughout the great majority 
of each of the two 2CK regimes. 

IV. QUANTUM PHASE TRANSITION 
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FIG. 11: T = spin-spin correlation functions (Si -82) (solid 
line) and (Si ■ S3) (dashed) vs J'/J, for fixed U — 7 and 
pJ = 0.1. Inset: magnification around the point of frustra- 
tion, where (Si • S2) = (Si ■ S3), occurring at J'/J = 1.032... 



respectively in the (— )-parity phase J' < to and — | 
respectively in the (-l-)-parity phase for J' > J^. That 
situation changes qualitatively on coupling to the leads, 
as illustrated in Fig. 11 where the T = spin correlators 
(Sl • S2) and (Sl • S3) are shown as a function of J'/J for 
the full model'*'''. While their behavior sufhciently deep 
in either of the 2CK regimes naturally accords with ex- 
pectations from the isolated TQD, they now vary contin- 
uously with J'/J: as seen clearly from the magnification 
of the crossover region in the inset. The value of J'/J 
for which ((S3 • S2) =) (Sl • S2) = (Si • S3) is the point 
of complete magnetic frustration, and as seen from the 
inset to Fig. 11 is naturally slightly renormalized from 
unity, to J^/J = 1.03216... in the example shown. 

We shall see in the following sections that the QPT 
occurs at Ea = J'^ — J' = 0, with a quantum critical 
point separating the two 2CK phases of distinct parity. 



When \E/\\ = \J — J'\ is sufficiently small however, 
neither TQD doublet alone dominates the low-energy 
physics, both must then be included in the low-energy 
trimcr manifold, and as shown in Sec. II B the low-energy 
behavior is no longer described by an effective single-spin 
2CK model, but rather by Eq. 13; with coupling between 
doublets of distinct parity embodied in the pseudospin 
raising/lowering terms. We show in the following that 
these additional terms in the effective Hamiltonian drive 
a quantum phase transition between the two 2CK phases, 
occurring at the point of inherent magnetic frustration. 

In the isolated TQD, frustration occurs at J' = J 
where the TQD doublets are degenerate, a trivial level- 
crossing 'transition' occurring in this case as = J is 
crossed. As a result, ground state properties in general 
change discontinuously across J^; exemplified eg by the 
spin correlation functions (Si • S2) and (Si • S3) which 
(as in Sec. II A) change abruptly at J'^ from — i and +j 



A. Physical picture of the transition 

Before showing NRG results, we give heuristic physical 
arguments for the behavior of the system in the vicinity of 
the QPT, based on the effective low-energy Hamiltonian 
^trans ^j^j^jj ^^^^j^ rpQp ^o^^let states are 

retained in the low-energy trimer manifold. H^ff^^ is 
cast in terms of spin-^ operators for both real spin (S) 

and pseudospin (7") for the local TQD Hilbert space, the 
distinct parity of the two doublcits being reflected in the 
= components of pseudospin. 

In Eq. 13 the exchange couplings ^{Ja^Jb) (= Jk±) 
are both > (AF), in which case it is readily argued (and 

confirmed explicitly by NRG) that the low-energy FP 
structure of H^g"^^ is independent of Jb- For simplicity 
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in the following we thus consider Js = in Eq. 13, ie 

^trans _ 1 J^g • (s^o + Srq) 

+ Jmix(t+ + f-)S ■ {SLO - SRo) (32) 

with £'a = J'c — J' (rather than the 'bare' = J — J' ^ 
allowing simply for the slight renormalization of away 
from the value = J in the isolated TQD, as above). 

For 7^ 0, a 'magnetic field' acts on the ^-component 
of pseudospin, and as mentioned in Sec. 11 B the low- 
energy physics is ultimately that of a single-spin 2CK 
FP (we return to it below). But now consider = 
in Eq. 32; ie J' = J'^, corresponding to the transition 
itself. In this case no field acts on the pseudospin, its 
'z-component' as such being arbitrary. Recalling that 
T+ + 7'~ = 27'^ in Eq. 32, and then performing a trivial 
rotation of the pseudospin axes from (a;, y, z) {z, x, y), 
Eq. 32 may be written as: 

Hlr^E^ = 0) = (i Ja + 2J^i.r^)S • s^o ^33^ 

+ {^Ja — 2JmixT^)S • Sro 

This Hamiltonian commutes with T^, so is strictly sep- 
arable into disjoint = +^ and — ^ sectors, ie the 
pseudospin is free. But in either given sector, the 
spin S clearly couples asymmetrically to the L and R 
channels, ie one has channel anisotropy. In the single 
spin-i 2CK model, channel anisotropy is of course well 
known4^'^^'5i'6i'*52 to destabiUse the 2CK FP: the spin-i 
is instead fully quenched by the channel to which it is 
most strongly AF exchange-coupled (and is decoupled 
from the second conduction channel/lead). The ulti- 
mate stable fixed point is then the single channel strong 
coupling (SC) FP, reached below a temperature scale 
T^^^ oc exp(-l//c)J>) with J> the larger of the L/R 
exchange couplings. 

In the present context, Eq. 33, the situation is then 
clear: in the = +5 sector the spin is wholly quenched 
by coupling to the L lead/channel (Jmix > 0), while 
for = — ^ it is quenched by coupling to the R lead; 
the temperature scale for quenching in either case being 
the one-channel Kondo scale T}f^ cx cxp(— 1//9J>) with 

The stable FP is clearly a doubled version of the SC 
FP, the 'doubling' refiecting the free pseudospin; with an 
associated T = residual entropy of Simp = ln2, and 
a finite uniform spin susceptibility Ximp(O) <^ l/Tjf'^ 
symptomatic of the 1-channel quenched spin-i. Since 
this FP is distinct from the stable 2CK FP arising for 
Ea ^ away from the QPT, it corresponds to a critical 
FP (CFP); with Eq. 33 the effective Hamiltonian for the 
quantum critical point (QCP) itself. 

The 'entire' T-dependence of thermodynamics at the 
QCP is also readily inferred, since the pseudospin is ubiq- 
uitously free and spin-quenching is characterised by the 
single, finite scale T^^. For T > one expects 



Simp{T) = ln4 (reflecting the free pseudospin and the 
free spin), with a free spin-i Curie law 2^Ximp(-^) ~ i' 
the 'SU{2) X SU{2y FP Hamiltonian here being Eq. 33 
with all exchange couplings set to zero. T ^ Tjf^ then 
sets the scale for spin quenching and the crossover to 
the CFP discussed above. Hence, at the QCP, the T- 
dependence of the entropy should be given by S'imp(T^) = 
hi2+Sl^^{T), with Sl^^{T) the entropy for a 1-channel 
Kondo model with Kondo scale Tjf^-. likewise the uni- 
form spin susceptibility should be Ximp(^) = xlmpi"^)- 

Having considered the QCP we return briefly to a non- 
vanishing I^^aI TfP^ , small compared to the ICK 
scale. In this case we expect flow to the CFP to be cut 
off at a characteristic 2CK scale, T ~ Tk, below which 
the system flows to the stable 2CK FP; and with Tk = 
Tk{Ea) vanishing as |£^a| — >■ 0, indicative of the QPT. 
We examine this via NRG in the following section, but 
note here that the considerations above would lead us 
to expect a crossover in SimpiT) from In 4 to In 2 on the 
scale T - T^<^^, followed by a crossover for T ^ Tk to 
the 5 In 2 characteristic of the stable 2CK FP. Likewise 
for Xinip(T') we anticipate the Curie law Ximp(^) l/4r 
for T > T^'^^, increasing to Ximp(^) ^ 1/T}f^ for 
T ^ Tjf^, before crossing over to the divergent 2CK 
behavior Ximp(T^) « HTk/T) for T < Tk- 

Finally, for a larger non-zero |£'a| ^ Tjf^, no RG 
flow in the vicinity of the CFP is expected, and the Tjf^ 
scale is then irrelevant: below T ~ \Ea\ ^ T]f^ one or 
other of the pseudospin ^-components is simply frozen 
out, >S'inip(^r) crossing over from In 4 to In 2 (here indica- 
tive of a frozen pseudospin but a free spin-i), ^nd then to 
the stable 2CK FP value i In 2 on the 2CK scale T ~ Tk- 
The physics here is simply that of a single spin-i 2CK 
model, as already considered in the preceding sections. 



B. Critical fixed point: NRG 

The physical arguments given above imply that the 
NRG energy levels associated with the critical FP itself, 
should consist of both a set of single- channel strong cou- 
pling FP levels, reflecting spin-quenching to one lead (as 
arises for an AF-coupled single-channel spin-i Kondo 
model); and a set of levels for a free conduction band, 
reflecting decoupling of the spin from the other lead (as 
arises at the trivial weak coupling FP for a ferromagnet- 
ically coupled single-channel spin-^ Kondo model^). 

Before proceeding we simply demonstrate that this is 
indeed the case. The left panel of Fig. 12 shows the low- 
est NRG energy levels of the full TQD model at the tran- 
sition (J' = J^), as a function of iteration number, N . 
The CFP levels are well converged by = 40 (flow to the 
CFP from the SU{2) x S'C/(2) FP beginning for A^ w 26). 
For comparison we also show the lowest NRG levels for 
a representative ICK model with AF coupling (middle 
panel) and with ferromagnetic coupling (right panel). 
The CFP level structure can indeed be seen to comprise 
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FIG. 12: NRG energy levels for even iteration number, N. The lowest 4 levels are shown in each charge/spin subspace for: the 
full TQD model precisely at the transition [left panel], a ICK model with pJx = +0.056 [middle] and pJk = —0.056 [right]. For 
the full model, subspaces shown are (Qi, Qr, S'') = (0, 1, 0)/(l, 0, 0) [solid lines], (0, 1, 1)/(1, 0, 1) [dashed], and (0, 3, 0)/(3, 0, 0) 
[dotted]. For the ICK models, the subspaces are (Q, S'^) = (0, 0) [solid lines], (1, i) [dashed], and (0, 1) [dotted]. 



both sets of strong coupling (middle panel) and weak 
coupling (right panel) levels. Indeed for the full TQD 
model, levels in the subspaces {Ql,Qr,S^) = {x,y,S^) 
and {y,x,S^) arc degenerate, so overall L/R symmetry 
is naturally unbroken, and the strong coupling and weak 
coupling levels form symmetrically in each channel. 

C. Thermodynamics of the transition 

We turn now to thermodynamics in the vicinity of the 
transition, obtained via NRG for the full TQD model; 
considering temperatures T <C F (ie without comment 
on non-universal high-temperatures, where the same be- 
havior as in Sec. Ill A naturally occurs). 

Fig. 13 shows the T/F-dependence of the entropy 
Simp{T). Results are given for fixed U = 7 and pJ = 0.1, 
varying J' as the transition occurring at J' = deter- 
mined as above from the spin correlators — is approached 
from either side, according to E/^ = (J^ — J') = ±ATg-; 
with results shown for A = 10" and n = 0, —1, —2, —3, —4 
[lines (a-e) in Fig. 13], chosen to approach progressively 
the transition, itself occurring at A = [shown as line (f)]. 
Here, T^/F = 10"^ is the value, at J' = J^, of the Kondo 
scale obtained from Txiy^piT) as in Sec. Ill A; and 
we will in fact shortly identify Tj^ as corresponding to 
the effective single-channel Kondo scale T]f^ discussed 
in Sec. IV A and associated with the QCP. Lines (a) in 
Fig. 13 are thus for J' being iT]^ away from the transi- 
tion; while (b-e) show the behavior for | — J'j <C as 
the QPT is approached. 

Sim.-p{T) self-evidently shows a QPT: one sees clearly 
a low-energy scale Tk = T^^ in the vicinity of the tran- 
sition, reflected in the crossover to the stable 2CK FP 
value S'imp = 5 In 2 (and determined in practice via 
5'imp(T'|) = 0.55 as in Sec. Ill A); with Tk vanishing 
precisely at the transition itself (line (f)). 

The obvious first question is how = Tk{E^) van- 
ishes as Ea = — J' — >■ on approaching the QPT. The 



answer is 

Tk '^'^^'■'^ AlJ'^-J'l" :iy = 2 (34) 

with exponent u = 2, and common amplitudes A on ap- 
proaching J'^ from either side. This is shown in the inset 
panel (c) to Fig. 14 (itself discussed further below). The 
solid lines therein show Tk/T vs \J'^ — J'\/T^ on a log-log 
scale, as the transition J' — >■ is approached from both 
sides. The dashed line has a slope v = 2, onto which the 
Tk fall clearly for | - J'j < T^; so that T^ is as such 
the 'boundary' scale below which Eq. 34 holds and the 
critical regime is entered. We add also that an exponent 
v = 2 reflects the non-trivial nature of the transition: 
for a first order level-crossing transition, as arises^'^ if the 
TQD is instead coupled to a single lead/channel at dot 
2 (see Fig. 1), one instead finds^'' v = \. 

Consider now the QCP itself, J' = J[. (line (f) in 
Fig. 13). From the arguments given in Sec. IV A we an- 
ticipate S'r^^^^/iT) = ln2 + Sl^^^{T), where the explicit 
In 2 reflects the free pseudospin and Sl^^{T) is the en- 
tropy for a sm^iZe-channel spin-^ Kondo model. This is 
indeed verified in Fig. 13, where In 2 + S'i^^p^(r) is calcu- 
lated explicitly for a 1-channel Kondo model with T}f^ 
chosen such that TjP^ = T^. It is seen (dotted line) to 
coincide perfectly with S^^^{T) for all T ^ T; includ- 
ing the In 4 entropy plateau reflecting the free pseudopsin 
and spin (the 'S'C/(2) x SU{2y FP of Sec. IV A) , and the 
crossover at T T]f^ = T^ to the CFP characterised 
by S'imp = ln2, where the pseudospin remains free but 
the spin is wholly quenched. 

Moving away from the QCP, consider now the sequence 
(b — >■ e) in Fig. 13, as the transition is progressively ap- 
proached. In all cases 5imp(T) follows the QCP S9^/{T) 
all the way from the In 4 regime, through the crossover at 
the common temperature T ~ T^ = Tjf^ (which is of 
course finite at the transition), to the In 2 symptomatic of 
the CFP; before ultimately descending to the stable 2CK 
FP with Simp = 5 ln2 on the scale T ~ T/f = T|. Since 
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FIG. 13: Entropy 5'imp(r) vs T/T on progressively approach- 
ing the QPT as: Ea = {Jc - J') = iATjf , with T^/T = IQ-^ 
(as discussed in text) and results shown for A = 10" with 
n = 0,-1,-2,-3,-4 (lines (a)-(e) respectively). The QCP 
itself (A = 0) is shown as line (f), while the dotted line shows 
In 2 + SlSp (T) calculated for a single-channel Kondo model 
with TjP^ = Tj^; the two coinciding for all T/F < 1. Inset: 
results in (b)-(e) rescaled (solid line) in terms of the vanishing 
Kondo scale Tk = T^, showing universal crossover from the 
CFP to the stable 2CK FP; and contrasted to that arising in 
the single-spin 2CK model (dashed line). 



Tff vanishes as the transition is approached, Si^p{T) 
should exhibit universality in terms of T /Tk- That it 
docs is shown in the inset to Fig. 13, where all lines (b- 
e) collapse to a common scaling curve (solid line). This 
scaling is of course characteristic of flow from the vicin- 
ity of the critical FP to the stable 2CK FP. As such it is 
thus expected to be distinct from the universal crossover 
in the single spin 2CK model (Sec. Ill A), from the In 2 
entropy plateau associated with the free spin-i local mo- 
ment FP to the stable 2CK FP; indeed confirmed in the 
inset to Fig. 13, where the latter is shown as a dashed 
line. Note also in this regard that lines (a) in Fig. 13, 

' ary 



±T]^, constitute in effect the bounch: 



for - J' 

between effective single-spin 2CK behavior (arising for 

i^aI = u: - j'l > T£ 



J' 



= j.^ j and the critical regime 
< T£-; as reflected in the direct flow of iS'iinp(T) 
from In 4 to | ln2 occurring at that point. 

Fig. 15 shows the T-dependence of the uniform spin 
susceptibility, T^Ximp(T) vs T/Tf^ (with T^/T = 10"'' 
as above), on approaching the transition according to 
Ea = J'c - J' = ATf ; but now with A = 10"/^ and 
n = 1,0,-1,-2,-3 (solid lines, in order of decreasing 
slope), showing as such the region where RG flow near 
the CFP is first observed (roughly between lines (a) and 
(b) in Fig. 13). The figure also shows the behavior at the 
QCP itself (A = 0, bottom solid line), as well as the spin 
susceptibility for a single-spin 2CK model (dashed line). 

(T) 



For the QCP itself, we expect X^Z^iT) = 



ICK 
Aimp 



from the arguments of Sec. IV A; with xlmp^ (T) the spin 
susceptibility for a single-channel spin-i Kondo model 
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FIG. 14: Variation of the Kondo scale Tk = T|: Tk /T vs 
(J' - J'c)/T^ (with fixed U = 7 and pj = 0.1, for which 
Tk/T = 10^^). Inset (b) shows an expanded view in the 
vicinity of the transition (solid line); inset (c) shows the same 
results (solid) on a log-log scale, together with the power- 
law decay Tk oc — J'p (dashed line) as the transition is 
approached. In both insets, is also shown (dotted lines), 
and remains finite at the transition. 



with Kondo scale Tj!f ^ = T^, such that at T = in par- 
ticular XimpC^ = 0) is characteristically finite. That is 
indeed verified in Fig. 15, where the result for the single- 
channel Kondo model is shown as a dotted line. 

For n = 1 by contrast, ie J'^ — J' « 2T^, the sus- 
ceptibility follows essentially perfectly the dashed line; 
thus being described by a single-spin 2CK model, pre- 
cisely as arises deep in either 2CK phase of the full TQD 
model (see Fig. 4), and characterised by the low-T log 
divergence T^Xim-p{T) oc ln(TR-/r). On moving closer 
to the transition (decreasing n), the susceptibilities pro- 
gressively 'fold on' to the QCP result over an increas- 



ingly wider T/T^-range; indeed even for — J' = O.IT^ 
(n = —4, not included in Fig. 15), the resultant suscep- 
tibility is barely distinguishable from the QCP line over 
most of the T/Tj^ range shown. As seen from the figure 
however, in all cases except the QCP itself, the ultimate 
low-T behavior is the log divergence expected for the sta- 
ble 2CK FP, TI^Ximp{T) ~ x\n{TK/T)] with an ampli- 
tude X visible from the gradients in Fig. 15 and seen to 
diminish steadily as the transition is approached. While 
numerical accuracy prevents a definitive determination of 
X much closer to the transition, extrapolation of the re- 
sults in Fig. 15 (as shown in the inset) suggest x vanishes 
as X oc Tk/T'}. (oc (1 - J'/J'^f from Eq. 34). 

Finally, as discussed only partially above (see Eq. 34), 
we return to the evolution of the Kondo scale shown in 
Fig. 14; the main panel of which shows Tk = T^ vs 
{J' — Jc)/T^ as determined from the T-dependcnce of 
the entropy (Fig. 13). Inset (b) to Fig. 14 shows an 
expanded view of Tk in the vicinity of the transition 
(I J' - J'c\/T^ < 2), together (dotted line) with the scale 
T^ determined in practice from the spin susceptibility 
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FIG. 15: Uniform spin susceptibility T^Ximp{T) v.s T/T^ for 
fixed U = 7 and pJ = 0.1, varying J' as: Ea = {J'c — J') = 
XT^ with A = 10"/* and n = 1,0,-1,-2,-3 (solid lines, 
in order of decreasing slope). The QCP itself (A = 0) is 
the bottom solid line. Also shown (dotted) is a ICK model 
chosen such that Tjf^ = TJ^, which describes the QCP; and 
a single-spin 2CK model with Tk = (dashed). Near the 
transition Ximp(T) diverges as Tk Ximp(r) ~ x \n{TK/T), with 
amplitude x ~ Tk/T^ as shown in the inset (see also text). 



via^-* T^Ximp(T^) = 0.07 (as employed in Sec. IIIA). 
While Tk = in inset (b) naturally shows the vanish- 
ing of the Konclo scale as the QPT is approached, 
is seen to remain finite at the transition, as seen also in 
Fig. 10 of Sec. inc. This is precisely as it should be. 
The 'true', ultimately vanishing Kondo scale Tk is of 
course evident in the evolution of both Simp{T) (Fig. 13) 
and XimpiT) itself (Fig. 15). But TxXi„p(T) — from 
which T^ has been defined — vanishes as T ^ (the 
explicit factor of T 'killing' the low-T log divergence in 
XimpCT") itself). Indeed, as readily inferred from Fig. 15 
(and verified by explicit calculation), T'Ximp(^) dof^t- to 
the transition (|Jc— J'l/T^ < 0.1 or so) is indistinguish- 
able from the QCP behavior Txj^^ (T); and as such is 
naturally characterised by the finite temperature scale 
Tjf^ = T^. The above, important distinction between 
Tk and is however applicable only close to the transi- 
tion. As seen clearly from inset (b) to Fig. 14, the scales 
Tk = T^ and coincide away from the transition (in 
practice for | — J'\/T^ ^ 2 or so); as noted and used 
in Sec. IIIA when considering the model deep in either 
of the 2CK phases. 



D. Dynamics of the transition 

Here we focus again on the single-particle spectrum 
Di{lo) of dot 1 (or equivalently, £>3(w) for dot 3). The 
T = spectrum wTDiluj) vs co/T is shown in Fig. 16 for 
U = 7 and pJ = 0.1; varying Ea = J'c^ J' on approach- 
ing the transition from the (-I-) parity phase J' > J'^, 



according to Ea = ^^T^ with A = 10"^^ and integral 
n = +4 — > —4 for lines a— >-i respectively. The spectrum 
at the transition itself (A = 0) is indistinguishable from 

case (i). 

Two particular spectral features should be noted. 
First, the spectrum at the Fermi level is pinned at 
wTDi{w = 0) = 5 in all cases, including the transition 
itself. Second, the width of the Kondo resonance clearly 
increases as the transition is approached, being of order 
T^ = Tlf^ at the transition (with T^/T = lO""* as 
usual); in which sense the vanishing Tk scale evident eg 
in the evolution of the entropy (Fig. 13) does not show 
up in single-particle dynamics. 

That -kTD\{lo = 0) should be | away from the tran- 
sition is (as in Sec. IIIB) a natural consequence of the 
stable 2CK FP that ultimately arises (Sec. IV C). To 
understand why ■kTD-i{uj = 0) = ^ at the transition it- 
self, first recall the general result from Sec. IIIB that 
7rrDi(a;) = — 7r/9ylmt^(a;), with ta{i-u) the t-matrix for 
scattering of electron states in the a = L,R lead. As 
discussed in Sec. IV A, the effective Hamiltonian for the 
QCP (Eq. 33) is strictly separable into disjoint = =^5 
sectors, the two sectors possessing common eigenval- 
ues. In this case it is straightforward to show that 
ta{uj) = l[ta{T' = +^;co') +t«(T^ = -^;w)], where 
ta{T^;oj) is the ia-matrix calculated for a fixed in 
the QCP Hamiltonian Eq. 33 (with fL(r^ = ±i;w) = 
tR{T^ = T^; w) such that tt(a;) = tB.{u)) overall, as sym- 
metry requires generally). Consider then tL(w). As ex- 
plained in Sec. IV A, at the critical fixed point the spin-^ 
S is quenched in a one-channel fashion, by AF coupling to 
the L-lead for = +51 but to the i?-lead for T^' = —\. 
In consequence, tL{T^ = -I-5; w = 0) is equivalently the t- 
matrix, t^^^{uj = 0), for a one-channel AF spin-^ Kondo 
model. But t^lT^ — — 5;^ = 0) = by contrast, since 
for = — i the spin is quenched by coupling to the 
i?-lead, while the L-lead is entirely decoupled from it 
(so no spin-scattering of electrons in the L-lead can oc- 
cur). Hence, tL{ui = 0) = \t^^^{uj = 0). For a one- 
channel spin-i Kondo model itself, —TTprphid}'~^^{uj) = 
ttT D^'^ ^ (uj) , with D^'^^ {(jj) the single-particle spectrum 
of a one-channel, single-level Anderson impurity model 
in the singly-occupied Kondo regime of the model (as is 
physically obvious, but follows more formally via argu- 
ments directly analogous to those given in Sec. IIIB). 
Overall we thus infer that 7tTDi{0) = — 7r/9ylmti(0) = 
^nVD^^^iO). But nTD^^^iO) = 1 by virtue of the 
Friedel sum rule^; whence 7rri3i(0) = | arises, as indeed 
found, see Fig. 16. 

While the arguments above apply to uj ~ (and as 
such the CFP), one naturally expects the frequency de- 
pendence of the Kondo resonance in nVDi (oj) at the QCP 
to be that of ^ttTD^^'^ (i^). That this is so is demon- 
strated in the inset to Fig. 16. The scaling spectrum 
^ttTD^^'^ (oj) for the Anderson model (as a function of 
ijj/Tjf^) is shown as a dashed line, and compared to 
7rrDi(w) vs uj/T1^ for the fuh TQD model at the QCP 
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FIG. 16: TrFDiiuj) vs cu/r, varying Ea = JL - J' as the 
transition is approached from the (+) parity 2CK phase: 
Ea = -AT^ with A = 10"^^ and n = +4^-4 (a^i respec- 
tively) . The QCP spectrum (A = 0) is indistinguishable from 
case (i). Inset: QCP scaling spectrum vs oj/T^ (solid line, in 
practice coinciding with all spectra for A < 10"'^). The QCP 
spectrum is seen to be the scaling spectrum ^-kT D^'-' (oj) 
(dashed line) for a single-channel Anderson impurity model. 
The dotted line shows the distinct scaling spectrum of a 
single-spin 2CK model (onto which cases (a)-(c) scale). 



FIG. 17: Local dynamic susceptibility T^Xii'^) ^/Tk, 
varying Ea = J'c ^ J' on approaching the transition from the 
(-j-parity 2CK phase: Ea = AT|- with A = 10"/^ and n = 
+4 —4. TkXiCO) vanishes as the transition is approached. 
The dashed line, T^xi('^) ^ '^i shows the low-w Fermi liquid 
behavior of the single-channel Anderson model arising at the 
QCP (A = 0). Inset: Xi{cu) = xi{^)/{MM{St; St))^=o]^) 
vs io/Tj^ at the transition (A = 0), showing recovery of the 
Korringa-Shiba relation Xi{uj = 0) = 1. 



(solid line); the two coincide perfectly. In particular, 

the low- frequency spectral behavior {\u>\/T^ ^ 1) at the 
transition is thus of quadratic Fermi liquid form^, 

nTD.iu^) ~ i[l-a'(u;/r^)2] (35) 

with a' a constant of order unity. This is of course in 
marked contrast to the behavior arising deep in either of 
the 2CK phases (Sec. lllBl). The scahng spectrum in 
that case is also shown in Fig. 16 (inset, dotted line), and 
instead exhibits the characteristic square-root frequency 
dependence of Eq. 27. 

Cases (a)-(c) in Fig. 16 are relatively deep in the (-|-)- 
parity 2CK phase (|£^a| > ^OT^). Their dynamics are 
thus in essence those of the single-spin 2CK model dis- 
cussed in Sec. IIIB: when scaled in terms of their Kondo 
scale Tk, the spectra all 'collapse' onto the 2CK scaling 
spectrum shown in the inset (dotted), departure from 
such occurring only by non- universal scales oj ~ \Ea\ 
(where a subsidiary peak arises, see main figure). Close 
to the transition by contrast exemplified by case (i) 
in Fig. 16 (|^a| = IQ-^T^) — the single-particle dy- 
namics are indistinguishable from that of the QCP, as is 
physically natural. Equally naturally, spectra (d)-(h) do 
not conform to either limiting form (QCP or 2CK), but 
instead represent crossover behavior between the two. 

Finally, we consider the T = local dynamic suscepti- 
bility for dot 1 in the vicinity of the transition. Fig. 17 
shows T^Xii^) ^f^K' varying = J'^ — J' on ap- 
proaching the transition from the (— )-parity 2CK phase 
J' < J'^, according to = AT^ with A = 10"/^ and 



n = -h4 H> -4. The cases n = 2, 3,4 (with ^a > lOT^) 
are quite deep in the 2CK phase (where the Tk scale 
is roughly constant, see Fig. 14); so they are essentially 
equivalent to the 2CK scaling curves of Fig. 9 (inset). 

As the transition is approached however, the height of 
the a; — > plateau in Fig. 17 is seen steadily to dimin- 
ish; the behavior found being of form Tj^Xii^ = 0) ~ 
Tk/T^, with the Kondo scale Tk 0. This is readily 
understood: as mentioned in Sec. IIIBl it is known^^'^^ 
that at the 2CK FP, Xi(0) plateaus at a constant, itself 
proportional to the slope of the low-temperature log di- 
vergence of Ximp(^)i ^^"i i'^ I^ig- 15 we showed the latter 
to vanish cx Tk 0. Hence Xi(0) = at the tran- 
sition itself. As expected from the nature of the QCP, 
the leading low-w dependence of Xi i'^) is then the Fermi 
liquid behavior characteristic of the single-channel An- 
derson model^'^^, Xii^) shown in Fig. 17^"^ (dashed 
line); and Xi(<^) is also seen to contain an absorption 
centered on w = (= T]f^), likewise characteristic of 
the Anderson model^'^^. 

The inset to Fig. 17 shows Xi{u]) = 
Xi(w)/(2w[Re((5f;.Sf))„=o]2) vs w/T^ at the tran- 
sition itself. For the single-channel Anderson model, 
the bj ^ behavior of Xi{ijj) is given exactly by the 
Korringa-Shiba relation^'^^ Xi{(jj = 0) = 1; seen from 
the figure to be well satisfied in practice (to within a few 
%), again confirming the physical picture of the CFP 
discussed in Sec. IV A. 
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SU(2) X SU(2) 



LM(-) 



FIG. 18: Evolution of the 2CK scale for the reduced model 
Eq. 36: Tk /D vs Ea/D (with D the conduction bandwidth), 
for fixed pJmix = 0.075. Inset: Tk/D vs \Ea\/T^ on loga- 
rithmic axes, with Tk ~ I-E^a]^ on approaching the transition. 



V. REDUCED MODEL FOR THE TRANSITION 

In the preceding sections, the effective low-energy 
model Eq. 13 (or Eq. 32) has been important in un- 
derstanding the behavior of the full TQD system in the 
vicinity of the QPT. We have also performed direct NRG 
calculations on the effective low-energy model itself, vary- 
ing independently the bare parameters Ja,Jb, Jmix and 
-Ba (or -Ea); and now comment briefly on the results of 
such. Specifically, we here consider explicitly the effec- 
tive model with J a = = Jb, which we have confirmed 
describes the same physics as arises with gcncrically non- 
zero Ja or Jb (reflecting the fact that the key element in 
the effective Hamiltonian Eq. 13 is the pseudospin rais- 
ing/lowering term, iJmix below, which in switching the 
parity of the TQD states effectively drives the QPT) . 

Wo thus consider the reduced model H^ed = Hq + 

-fornix + -ffmag) with Hq 

the Hamiltonian for the leads and 



J ■ (T^ 

mix — <^mixV ' 



mag 



T )S • (slo - sro), (36a) 

(36b) 



The underlying FPs of i?red are readily inferred. The 
SU{2) X SU{2) FP corresponds to Jmix = = £'a in 
Eq. 36, and as such to both a free spin and a free pseu- 
dospin; this of course is the high-temperature FP of i?red, 
generating eg the associated In 4 entropy. Two LM FPs 
also arise, corresponding formally to Jmix = in Eq. 36 
(and hence a free spin), and E/\ — >■ ±oo. For E/\ — >■ -|-oo, 
the = pseudospin component {ie the (— )-parity 
TQD doublet) is frozen out, and we designate this as the 
'LM(-)' FP; while for Ea -oo the T"" = +^ compo- 
nent is frozen out, corresponding to a LM(-|-) FP. 

The critical FP arises for Ea = 0. In this case (see 
Sec. IV A) , a trivial rotation of the pseudospin axes gives 



FIG. 19: Schematic illustration of RG flow between the FPs 
(circles) of the reduced model. See text for discussion. 



the QCP Hamiltonian H^'~^^ 



'red -Ho+ 2J„iixr^S • (SiO - 

Sflo)j with the quantum number conserved and the 
Hilbert space thus separable into disjoint = ±i sec- 
tors. As discussed in Sec. IV A, the CFP itself corre- 
sponds to the spin being quenched in a one-channel fash- 
ion by coupling to the L lead in the = +^ sector, 
and to the R lead for = — i. More formally, on writ- 
ing if^^f = J?o + 2P-S ■ (J,ti,Sio - J,'!;jBn) with 
■^mix = "^mix = -^mix > in the 'bare' H^^"" — the CFP 
corresponds to J^^^ — )• oo and J^j^ = in the — +^ 
sector, with J^;^ = and J^;^ — )• oo for = — ^. 

Finally, for E^ 7^ 7^ ^mix, the higher lying compo- 
nent of the pseudospin is naturally frozen out on suffi- 
ciently low energy scales T <C |-E'a|- Virtual excitations 
to the upper pseudospin sector still arise, however; and 
treating the mixing term Eq. 36a perturbatively (strictly 
valid for large |£;a|) gives H^^^ = Ho + H^%, with 



rreff 
■"hyb 



( ■^mix \ 
2EaJ 



S • (slo + Sflo). 



(37) 



This is a model of 2CK form in the reduced Hilbert space 

of the lowest pseudospin component ((+)-parity for E^ < 
and (— )-parity for Ea > 0), for which the ultimate 
stable FP is of course the infrared 2CK FP. 

The overall low-energy FP structiire of the reduced 
model is thus indeed that of the full TQD model (we 
discuss its RG flow diagram below). On studying Eq. 36 
directly with NRG the transition occurs as expected at 
Ea = 0, and the above FPs are indeed observed eg in 
the T-depcndence of thermodynamics. 

Fig. 18 shows the resultant 2CK scale Tk/D vs Ea/D 
for systems with common pJmix = 0.075. Tk is deter- 
mined in the usual way from the entropy, and foimd to 
be non-zero for any Ea ^ (Fig. 18); and we identify 
Tf. (oc exp(— l/pJmix)) as the crossover scale from the 
SU{2) X SU{2) FP to the CFP. The inset to Fig. 18 shows 
Tk/D vs \Ea\/T^ on logarithmic axes. Just like Fig. 14 
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for the full TQD model, the point Ea = can be seen 
to separate two distinct regimes. For |£'a| ^ T^i ^Y^' 
tems may be considered 'deep' in the 2CK phase (which 
situation constitutes the majority of the main panel of 
Fig. 18, and indeed the majority of the parameter space 
of Eq. 36). For ^ by contrast, systems flow close 
to the CFP, and are 'near' the transition. The vanish- 
ing Tk scale associated with the transition is found to 
be characterized by a power-law decay Tk ^ I-BaT with 
u = 2,as found in the full TQD model (Sec. IV C). We 
also add that since the reduced model Eq. 36 lacks a di- 
rect coupling of form S • (s^o + Sj^q) correlated to the 
pseudospin (c/ Eq. 13), the variation of Tk is symmetric 
in Ea; as evident in Fig. 18. 

A schematic illustration of the RG flow is given in 
Fig. 19, in the (.Ea , Jmix)-plane. The FPs are indi- 
cated by circles, the arrows representing RG flow be- 
tween them. At the highest temperatures, the system 
is always described by a free spin and a free pseudospin 
(the SU{2) X SU{2) FP, J^ix = £^a = in Eq. 36). 

Considering now the case of J^ix = 0, the spin al- 
ways remains free; but the upper pseudospin component 
is frozen out for any E^ 0, whence as T — )■ the sys- 
tem is described by either a (-I-) or (— ) parity LM FP. By 
contrast, for E^ = 0, any Jmix 7^ drives the system to 
the CFP as T — >■ 0. In the general case where \Ea\ 7^ 
and Jmix 7^ 0, the system always flows ultimately to a 
2CK FP (the parity of which depends on which compo- 
nent of the pseudospin lies lowest, ie whether Ea > or 
< 0). In this case RG flow first approaches the LM FPs 
if |i?A| > T^, or the CFP for \Ea\<S:T^, before finally 
flowing to the appropriate stable 2CK FP. 

And as evident from the results of Sec. IV {eg Fig. 13), 
the RG flows depicted schematically in Fig. 19 are just 
those arising at low-energies in the full TQD model, in 
the vicinity of, and at, the quantum phase transition. 



VI. CONCLUSIONS 

We have studied in this paper what is arguably the sim- 
plest 2-channcl model in which to understand the conse- 
quences of local frustration arising from internal degrees 
of freedom: a triple quantum dot ring, with dots mutu- 
ally coupled by antiferromagnetic exchange interactions, 
and tunnel-coupled symmetrically to two metallic leads. 
While important aspects of the general model have been 
considered before^^'^^, the underlying physics is found to 
be even richer than hitherto uncovered. 

Two distinct 2CK phases arise due to the mirror sym- 
metry in the problem, each displaying classic non-Fermi 
liquid properties of the 2CK flxed point below a charac- 
teristic 2-channel Kondo scale. But although 2-channel 
Kondo physics predominates in the underlying parame- 
ter space, the parity-distinct nature of the 2CK phases 
means that a quantum phase transition between them 
occurs. Driven by varying the interdot exchange cou- 



plings, occurring at the point of inherent magnetic frus- 
tration, and characterised by a nontrivial quantum crit- 
ical point which we have explicitly identified and anal- 
ysed, the transition provides a striking example of the 
subtle interplay between internal spin and orbital degrees 
of freedom. 



Acknowledgments 

Helpful discussions with M. Galpin, C. Wright and T. 
Jarrold are gratefully acknowledged. We thank EPSRC 
(UK) for financial support, under grant EP/D050952/1. 

Appendix 

Here we specify the parity operator Pij which ex- 
changes the labels i and j for an arbitrary pair of orbitals 
First we transform canonically to even (e) and odd 
(o) orbitals, 

cl=^^{cl + cl) (A.la) 
4. =73(4 -4)- (A.lb) 

Exchanging i ^ j clearly has no effect on cl^, but 
cl^ -B- —cl^. With ho = hof+fioi the odd-orbital number 
operator we thus define: 

Pij = 2{no (A.2) 

Pij is self-adjoint and involutory {Plj = 1), and hence 
has eigenvalues ±1 only (specifically Pij = +1 for no = 
0,2 and Pij = -1 for rio = 1). From Eqs. A.l and A.2, 
it follows that 

[hj,cl]=0, (A.3a) 
{4„cL} = 0, (A.3b) 

and hence from Eqs. A.l and A. 3 

hA. = <^].Ko (A.4a) 

KA = 4Aj , (A.4b) 

showing that P^j indeed permutes the i, j labels. Using 
Eq. A.l in Eq. A.2 it is straightforward (if lengthy) to 
show that 

Pij = 2(ii • i,- - Si • S,) + \ni^{ni-a-l) 

«={«>j}><^ 

+ E(4s-a + H-c-)(2 -hi-hj) + \ 

(A.5) 

where li is an isospin operator (with components If = 
i(n. - 1) and /+ = 4^,^ and J" = (J+)t). 
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If a Hamiltonian is invariant to the permutation i ^ j, 
then [H, Pij] = 0, parity is conserved, and all states can 
thus be classified according to it. The isolated trimer 
Hamiltonian iJtri (Eq- 1) is clearly invariant under inter- 
change of the 1^3 dot levels, Pi^^. The fuU lead-coupled 
Hamiltonian is by contrast invariant to simultaneous in- 
terchange of 1 and 3 labels and left and right leads, the 



larger ('overall L o i?') symmetry of H implying that 
the relevant involutory permutation operator is 

PL,R = Pl,3l[PLk,Rk , (A.6) 
k 

such that [H, Pl^r] = 0. 
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